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We investigate theoretically the behavior of a one-dimensional interacting Bose-Fermi mixture 
with equal masses and equal repulsive interactions between atoms at finite temperature in the 
scheme of thermodynamic Bethe Ansatz. Combining the Yang- Yang thermodynamic formalism with 
local density approximation in a harmonic trap, we calculate the density distribution of bosons and 
fermions numerically by treating the radially and axially excited states as discrete and continuous 
ones, respectively. Our result from exactly solvable solutions may be used as a touchstone of ID 
interacting Bose-Fermi mixture for experimental data fitting where mean-field theoretical approaches 
fail. 
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I. INTRODUCTION 

Recent years have witnessed great development of laser 
cooling and optical trapping technology, with a remark- 
able achievement being the quasi one-dimensional (ID) 
quantum gases realized by tightly confining the particle 
motion in the other two directions to zero point oscilla- 
tion [l|, 12, S S H [ll • Meanwhile the Feshbach resonance 
may be used to manipulate the interparticle scattering 
length by simply tuning an external magnetic field, which 
enables us to explore the cold atomic gases in the full 
interaction regime from weakly to strongly interacting 
limit. In particular, a quasi-lD quantum gas of strongly 
interacting bosons has been observed in the so-called 
Tonks-Girardeau (TG) regime On the other hand, 

mixture of quantum degenerate gases form novel quan- 
tum many-body systems with rich phase structures. Par- 
ticularly interesting systems are the Bose-Fermi mixtures 
0, nil, which rarely occur in nature and have become 
experimentally accessible with the development of sym- 
pathetic cooling [1, . Theoretical investigations on the 
quasi- ID Bose-Fermi mixture have focused on their phase 
diagrams and ground-state properties in the scheme of 
Luttinger liquid theory [l^, [HI and Bethe ansatz method 

The ground state energy of a ID system of bosons with 
repulsive delta- function interaction was first calculated 
by Lieb and Lininger • Yang and Yang extended this 
Bethe ansatz method to finite temperature 40 years ago 
p^ . The Yang- Yang thermodynamic formalism, also 
known as the thermodynamic Bethe ansatz (TBA), al- 
lows to evaluate the thermodynamic properties of the 
ID system. It has triggered numerous further investiga- 
tions and has been generalized later to spin-1/2 fermions 
p^ . mixture of spin-1/2 fermions and bosons [l^, and 
two-component bosons (T6| . etc. The first direct com- 
parison between experiments and theory based on the 
Yang- Yang exact solutions was carried out last year [23| 



in the weakly interacting Bose gas on an atom chip and 
for a wide parameter range where conventional models 
fail to quantitatively describe in situ measured spatial 
density profiles. In view of the experimental progress, it 
is thus quite desirable to study theoretically the thermo- 
dynamic properties of the quasi-lD Bose-Fermi mixture, 
which might provide theoretical guidance on the poten- 
tial experimental implementation. From pure theoretical 
point of view, it is also very interesting to study how 
robust of the Bose-Fermi phase separation, which was 
predicted to appear in the limit of zero temperature p7| , 
against the temperature effect. 

In this paper, we combine the TBA and local den- 
sity approximation (LDA) to study the density distri- 
bution of a mixture of bosons and polarized fermions 
trapped in a harmonic trap at finite temperature. Here 
we have the temperature, particle number of bosons and 
fermions, and the interaction strength as variables. We 
find that the observation of Bose-Fermi phase separation 
requires even lower temperature attainable by present 
cooling techniques. 

The article is organized as follows. In Section II, we 
rederive the thermodynamic Bethe ansatz solution by 
means of Yang- Yang thermodynamic formalism for Bose- 
Fermi mixture where both bosons and fermions are spin 
polarized. Section III describes our numerical procedure 
on how to numerically evaluate the density distribution 
of bosons and fermions in a harmonic trap at finite tem- 
perature. We treat in a different way for radially and ax- 
ially excited states. Finally in Section IV, we analyze the 
low temperature behavior of the ID Bose-Fermi mixture 
for realistic experimental situation and make concluding 
remarks in Section V. 



II. THERMODYNAMIC BETHE ANSATZ OF 
BOSE-FERMI MIXTURE 
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We study a ID interacting Bose-Fermi mixture on a 
line of length L with periodic boundary conditions, de- 
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scribed by the Hamiltonian 



t,T,t, 



(1) 



where VPf,, are field operators for a boson of mass rrih 
and for a fermion of mass m/, and gt;,, g?,/ are boson- 
boson and boson-fermfon interaction strengths, respec- 
tively. Fermions are spin-polarized so that Pauli principle 
excludes their s-wave interaction [g^f ~ 0). This model 
is exactly solvable for equal masses and equal repulsive 
boson-boson and boson-fermion interaction strengths, 



rrih = rfif 



gbb = gbf = 9- 



(2) 



Although an exact solution is available only under condi- 
tions ([2]), deviations slightly from this integrable line are 
expected not to dramatically change the characteristic 
properties of the system, e.g., the phase separation. Fol- 
lowing Lai- Yang's original convention, we assume 2m = 1 
and h = 1 and write the Hamiltonian ([1]) in its first quan- 
tization form 



H 



N 

i—i ^ 



2c^ ^ {xi - Xj) , 



(3) 



with c — mg/fi^. Among the N particles there are 
M bosons and the rest of them are fermions. The 
many-body wave function is supposed to be symmet- 
ric under odd permutations with respect to indices i = 
{!,..., M} (bosons) and antisymmetric with respect to 
i = {M + 1, A''} (fermions). 

For periodic boundary conditions, Imambekov and 
Demler obtained the following set of Bethe Ansatz equa- 
tions (BAE) [13 



M 



exp {ikjL) = JJ^ 



13=1 



kj - Kp + ic/2 
Kp - ic/2 ' 



,N, (4) 



1 = 



N 

n 



ki — Aq, -|- ic/2 
ki — Aa — ic/2 ' 



1,...,M, 



(5) 



where the momenta fci, ...fcjv are a set of unequal num- 
bers, and spectral parameters Ai, ...Aa/ are the analogs 
of the momenta. It has also been proved that all solu- 
tions of ([3][5]) are always real, which simplifies greatly the 
analysis at finite temperature. 

For thermodynamics at finite temperature we use the 
BAE to derive a set of nonlinear integral equations, i.e. 
TBA equations, which describes the thermodynamics of 
the model at finite temperature. Taking logarithm of 
the BAE (jUlS]), we arrive at the following discrete Bethe 



ansatz equations 

kjL ^ 2t:Ij + {2k j - 2Ap) , 



M 



/3=1 



N 



271 Ja = {2h - 2Ka) , 



(6) 



with 9 {k) = — 2 arctan (fc/c). Here Ij and Ja are inte- 
ger or half integer quantum numbers (depending on the 
parity of M and N) , which play the role of quantum num- 
bers for the momentum k and spectral parameter A re- 
spectively. For a particular configuration, if an arbitrary 
quantum number is chosen, it is either occupied (in the 
set of quantum numbers for the system), called a root, 
or not occupied, called a hole. In the thermodynamic 
limit, the distributions of momentum and spectral pa- 
rameter become dense, and it is convenient to introduce 
the density functions of roots and holes, respectively. We 
denote with p (k) and ph (k) the density functions of the 
momentum k and holes, and with a (A) and ah (A) the 
density functions of spectral parameter A and its holes. 
They are defined by 

Lip{k)+phik))dk^dl, 
L {a {A) + ah (A)) dA ^ dJ. 

Differentiate equation ^ with respect to k and A sepa- 
rately, we obtain a set of coupled integral equations 

p{k) + ph{k) = ^ + ^ p K{k,K)a{A)dA, 
a(A) + ah(A) = ^ f K {A, k) p (k) dk, (7) 



where 



K{x,y) 



4c 



-\- A{x — y) 



The total number of particles and that of bosons per 
unit length can be obtained by integrating the density 
functions of momentum and spectral parameter as follows 

/OO f^OO 
p{k)dk, M/L= / a{A)dA, (8) 
-OO -/ — OO 

while the energy of the system per unit length is given 

by 



/OO 
k^p{k)dk. 
-OO 



With the help of the approach first introduced by Yang 
and Yang, the entropy of the present model at finite tem- 
perature is 



S/L = J [{p + ph) In {p + ph) - pin p - ph h\ph] dk 
+ I [{a + ah) \n{a + ah) - a\na - ah In ah] dA. 
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The Gibbs free energy of the model is then defined by 
F = E - TS - fiF {N - M) - hbM, where fip and hb 
are two Lagrange multipliers, and T is the temperature. 
In order to arrive the thermal equilibrium, we minimize 
the free energy with respect to the density functions p (fc) 
and cr (A) subject to the constraint ([7]). In addition, the 
numbers of fermions and bosons are kept to be constants 
respectively. It can be proved rigorously that fip and 
fiB are the chemical potentials of fermions and bosons, 
respectively. 

Applying the minimum condition SF = gives rise 
to the following nonlinear integral equations, i.e. TBA 
equations 

T 

e (fc) = -^F + k'^ - 

Ztt 

/oc 
K (fc, A) In (1 + exp {-ip (A) /T)) dA, 
-OO 

T 

{K) ^ llF - llB - TT 

/OO 
(A, /c) In (1 + cxp (-£ (A) /T)) dfc, (9) 
-OO 

where we have defined 

exp(e(fc)/T) = ph{k)/p(k), 
exp(^(A)/r) = a„(A)/a(A), (10) 

and the set of equations (O becomes 

/OO 
K{k,A)a{A) dA, 
-OO 

/DO 
K{A,k)p{k) dk. 
-OO 

(11) 

The density functions p (fc) and a (A) can be obtained 
by solving the above coupled integral equations (O and 
pT|) . The TBA approach described here is universal for 
discussing the thermodynamics of ID integrable model. 
Once T,c,p,F and pb are determined, all thermodynamic 
properties are known. For instance, the pressure and the 
free energy are 

F = -PL + MpB + iN ~ M)pF- (12) 

III. LOCAL DENSITY APPROXIMATION AND 
NUMERICAL PROCEDURE 

Starting from the solution to the TBA equations de- 
rived above for a uniform ID Bose- Fermi mixture in ther- 
mal equilibrium, in this section we aim to find the nu- 
merical results for the finite-temperature density distri- 
butions of bosons ub (x) and fermions n p (x) under the 
local density approximation (LDA). 



Firstly we describe our iteration process in solving the 
TBA integral equations. Due to the fact that the TBA 
equations are a set of coupled nonlinear integral equa- 
tions there exist no closed analytical solutions for them. 
Nonetheless, numerically this is in principle a well con- 
trollable problem (it becomes, however, quite complex 
for an increasing number of particle types) and we here 
solve the equations by iteration. The convergence of this 
procedure to a solution and its very existence have been 
investigated most naturally by means of the Banach fixed 
point theorem [2l[. 

The iteration process is as follows. For given T,c,pF 
and pb, we initialize e (fc) and ip (A) on the rhs of equa- 
tions ([9]) with the corresponding zero-temperature trial 
functions e*^"^ (fc) = —jiF + k"^ and (/S^^^ (A) = pp — Pb, t^g- 
spcctively. In a first step, we obtain e^^^ (fc) and (p'-^-' (A) 
on the Ihs of equations ^ and let them be the new 
trial functions. The scheme continues with updates to 
e("' (fc) and .^("^ (A) with n ^ 1,2,.... When the rel- 
ative error between e^"^ (fc) and e("+i) (fc) and that be- 
tween (y?^") (A) and (^("+^) (A) reach a small quantity, e.g., 
10^^", sufficient convergence is obtained and e'-"-' (fc) and 
(^("^ (A) are considered as the solutions of TBA equa- 
tions. We then put these solutions into equations (fTTj) . 
and meanwhile the initial trial density functions p (k) 
and cr(A) is set as p'-^^k) = l/27r (1 exp (e (fc) /T)) 
and (T^°) (A) = 0, respectively. With the same iteration 
process, we can obtain the solution of equations (|lip . 
the integration of which gives the particle density. Note 
that, to insure the accurateness of the integration in the 
equations Q and pTjl for very small interaction strength 
c, the integrand should be divided into more parts with 
interpolation method ffi'stly and then integrated numer- 
ically. 

Experimentally ID gases are usually achieved by trap- 
ping the atoms in a tight harmonic trap with strong 
transverse confinement and weak confinement along the 
axis, L0± 3> w//. The main parameters from the experi- 
ment [lOl on Yang- Yang thermodynamics of *^Rb atoms 
in the \F — 2,mp — 2} hyperfine state are adopted here. 
The Bose-Fermi mixture is trapped in a harmonic poten- 
tial with u>±/2'K = 3280Hz, w///27r = 8.5Hz. 

The numerical procedure is as follows: The total lin- 
ear density relies on two contributions, i.e. that from 
the radial ground state and that from the radially ex- 
cited states. The density of atoms populated in the radial 
ground state is obtained by solving the TBA equation, 
while the atoms in the radially excited states are treated 
discretely as an independent ideal ID Bose(Fermi) gas 
in thermal equilibrium. This is because our tempera- 
ture here is on the order of the radial level splitting, 
hu!±/kB = 157. AnK, so that the fraction of the atoms 
in radially excited states should not be neglected. Along 
the axis, the energy level splitting, hLO///kB = O.inK 
is so small that we can use the LDA to account for the 
axial potential via a continuously varying chemical po- 
tential p (x) = p — V{x). In this way, we obtain the total 
linear densities of bosons UBix) and fermions np{x) in 
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the magnetic trap which may be used to fit the experi- 
mental data from absorption images. 

For the radial ground state, the TBA equations for 
uniform gas can be applied locally to the trapped gas if 
the condition for the LDA arc met. One assumes that 
in slowly varying external harmonic trap the chemical 
potentials of bosons and fermions are changed into 



Ms (X) 



(13) 



Here, fiB and /xf are chemical potentials for bosons and 
fermions in the center of the harmonic trap. Replac- 
ing the chemical potentials in eqs. ^ by their LDA 
values, we can obtain numerically the density functions 
p{^iB, ^J'F,x,k) and a {^j,b, IJ-FjX^A) by iteratively solving 
the TBA equations Q together with the constraint pT|) . 
The integration of these density functions yields the axial 
density distributions 



TBA 



(x) = / cr(jUB,yUF,a;, A)dA 



(14) 



for bosons and 



^TBA 



(x) 



p{fiB, fJ-FTX^k) dk ~ rig^^ (x) (15) 



for fermions, which are the LDA revisions to their uni- 
form counterparts, i.e. eqs. ([5|). 

Similar strategies are applied to acquire the densities 
of bosons and fermions at the radially excited states. For 
bosons it can be expected that the interaction will signif- 
icantly affect only the distribution in the radial ground 
state, while the population in the radially excited states 
can be well described by the distribution of ideal Bose 
gas. Fermions in the ground state do not interact with 
each other due to the Pauli exclusive principle. Thus for 
fermions, the chemical potential fiF can be even larger 
than hujj_, with the population in the radially excited 
states described by the distribution of ideal Fermi gas. 
The radially excited states are (j + l)-fold degenerate, 
i.e. for each radial quantum number j > 1 there are j + 1 
excited states sharing the same energy. We treat each 
excited state as an independent ideal ID Bose or Fermi 
gas in thermal equilibrium with the chemical potential of 
the gas in the radial ground state given by 



f^B i^) 
f^F i^) 



Mb (x) - jhuj±, 
Hf (x) - jhuj±, 



(16) 



respectively. The density distributions in radially excited 
state (j) are 



(x) 



1 




1 






2"" exp 
1 


y 2m 


- Ms {x)^j 

1 


IkBT 


- 1 


'^'"^ exp 


y 2m 


- Mf (^)^ 


/kBT 


+ 1 



dk, 



dk. 



(17) 
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FIG. 1: (Color online) Bose and Fermi density distributions 
at high (a) and low (b) temperatures for weakly interacting 
mixture. The interaction strength is chosen as co = 0.30/im. 
Yang- Yang thermodynamics is clearly seen for lower temper- 
ature as a narrow density peak. Also shown are the ideal 
Bose-gas profile and that from a mean field theory in the 
Thomas-Fermi approximation with the same peak density. 



The total linear densities are given by summing over 
the TBA results for radial ground state p4l [T5| and the 
ideal-gas results for radially excited states (fT7|) 

oo 

UB {x) - nl^"^ (x) + ^ (j + 1) n^B i^) > 

OO 

nF{x) ^nl^^{x) + Y.^j + \)n'p{x). (18) 

3 = 1 



IV. LOW TEMPERATURE BEHAVIOR AND 
EXPERIMENTAL CONSIDERATIONS 

Here, as an example which is experimentally accessi- 
ble, we show the linear density of atomic clouds with 3D 
scattering length a equal to lOOas {as is Bohr radius) 
and the mass of the atoms, both boson and fermion, 
being chosen as that of ^^Rb. ID Bose-Fermi mixture 
have not attracted much attention until recently, when it 
became possible to realize such systems in experiments 
with cold atoms [l^]. Experimentalists care more about 
the heteronuclear Bose-Fermi mixture, such as ^^Rb-^"K, 
in which case the heteronuclear interactions can be tai- 
lored by means of Feshbach resonances. Bose-Fermi mix- 
ture system, however, may be composed of isotopes of 
atoms whose mass difference is very small, such as ^Li- 
^Li, or ®^Rb-*^Rb, etc. The exactly solvable case con- 
sidered here is relevant to current experiments, and can 
be used as a benchmark to check the validity of different 
approaches. 

We first consider the weakly interacting mixture with 
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Hj,/k^=0.075nK 
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T=0.14nK 




N =100,N =100 





-Boson,T=15.74nK 
Fermion,T-15.74nK 
Boson,T=1.574nK 
Fermion,T=1 ■574nK 
-Boson,T-0.5nK 
Fermion,T-0-5nK 
Boson,T-OnK 
Fermion,T=OnK 





x(nm) 




FIG. 2: (Color online) Bose and Fermi density distributions 
with fixed chemical potential in the trap center and for differ- 
ent interaction strengths. When the interaction strength in- 
creases, less number of bosons are needed in the gas to achieve 
the chemical potential, while the number of fermions in the 
mixture almost remain unchanged. 



the effective ID coupling strength expressed through 
the 3D scattering length a as c — 2niLLj±a/h if a 

We denote as cq the coupling strength 
for bare ^^Rb background scattering which is approxi- 
mately O.SO/iTO^^. In Figs. 1(a) and (b), we show the 
linear density of atomic clouds for Nb bosons (black solid 
line) and Np fermions (blue dashed line) in the mag- 
netic trap for different temperatures. Experimentally 
these data may be obtained by absorption imaging and 
integrating the atom number along z-axis. At high tem- 
perature (T ~ 0.44/iif ), population in excited states con- 
tribute a lot to the density and the result from Yang- Yang 
formalism for ground state is only a small fraction, hence 
is not visible. The Yang- Yang thermodynamics is clearly 
seen for lower temperature (T = O.lAfiK) as a narrow 
density peak for bosons, where both the ideal-gas and 
the mean field distributions in the Thomas- Fermi approx- 
imation fail to quantitatively describe the spatial density 
profiles. We also find that in the weakly interacting limit 
the number of fermions does not affect the density profile 
of the bosons very much. For instance, one can hardly 
discern the difference in the distributions of 5 x 10'^ bosons 
when we include in the mixture 5 x 10^, 5 x 10^ or 5 x 10"* 
fermions. This can be attributed to the relatively small 
number of fermions in the ground state. 

The interaction in the mixture, on the other hand, may 
be enhancedgreatly by the so-called confinement induced 
resonances [5| . Additionally, all interactions can be tuned 
using available Bose- Fermi Feshbach resonances [22l|. In 
Fig. 2 we illustrate how the bosonic density distribution 
will change with the increase of interaction strength. For 
fixed chemical potentials of Bose and Fermi gas, as in- 



FIG. 3: (Color online) Bose and Fermi density distributions 
at different temperatures for large interaction strength c = 
100 fim~ ^ . Phase separation appears at very low temperature. 



teractions get stronger, the density of Bose gas decreases 
very quickly, while the density of fermi gas keeps un- 
changed. This can be explained as follows. In determin- 
ing the density of the atomic gas, Pauli exclusive prin- 
ciple plays a more important role than the Bose-Fcrmi 
interaction gb / for atoms in the ground state (even in the 
strongly interacting limit) as a result of most fermions 
occupying the excited states. 

Imambekov and Demler [l3] predicted the existence of 
the Bose-Fermi phase separation at very low tempera- 
ture and very strongly interacting limit, i.e. the relative 
distribution of bosons and fermions changes with interac- 
tion and the Fermi density shows strong nonmonotonous 
behavior for strong interactions. They firstly get the 
magnon energy spectrum for large 7( 7 — c/n and n 
is density) and then use the local density approximation 
and energy spectrum to obtain the density distribution. 
There is, however, no obvious signature of this phase sep- 
aration in our Fig. 2 which is plotted for T = OAA^K 
even for very large interaction strength. Clearly the 
temperature smears off this many-body quantum effect. 
To observe the phase separation, we need further cool 
the atomic gas to even lower temperature. In Fig. 3, 
we compare the density profiles for a mixture of 100 
bosons and 100 fermions in the strongly interacting limit 
c = 100/im~^ for different temperatures, some of which 
are beyond the current experimental reach. One clearly 
sees that the phase separation appears at nK tempera- 
tures. Also shown are the zero temperature result which 
is obtained analytically from Bethe ansatz method p7| . 

Finally, let us see how the number of fermions would 
affect the Bose density distribution. In the weakly inter- 
acting limit, we do not observe significant change in the 
density profile of the bosons when we change the number 
of fermions. However, when both the strongly interacting 
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FIG. 4: (Color online) Modification of bosonic density by dif- 
ferent number of fermions in the mixture for large interaction 
strength c = 100 fj.m~ ^ . In ultra-low temperature (nK) few 
bosons would populate on the radially excited states. For 
more fermions in the mixture, bosons in the trap center move 
to the edge. At slightly higher temperature, the total den- 
sity remains unchanged. Here "g" stands for density profile 
of atoms in the ground state. 

and ultracold conditions are met, the presence of more 
fermions drastically alter the density of bosons. We show 
in Fig. 4(a) the modification of bosonic density by differ- 
ent number of fermions. More fermions produce stronger 
Fermi pressure, which tends to flatten the Bose density 
when the interactions between bosons and fermions are 
strong. We also noticed that the ground state density 
(dashed lines) coincides almost perfectly with the total 
density (solid lines) for Np = 10 — 100, indicating that all 
bosons Nb — 100 occupy the radial ground state. One 
can see the distinction appears for even larger number of 
fermions, e.g., for Np = 200, the fattest two curves do not 
match each other any more. Temperature again blurs the 
quantum effect from increasing the number of fermions. 



i.e. the total density remains unchanged although that 
of atoms in the ground state indeed decrease for more 
involved fermions as can be seen in Fig. 4(b). 



V. CONCLUSION 

In summary, we have theoretically studied the quasi- 
ID system of Bose-Fermi mixture with equal repulsion 
between atoms at finite temperature. For the mixture 
system in a harmonic trap, we calculate the density dis- 
tributions of bosons and fermions numerically by using 
the combination of LDA and TBA and treating the radi- 
ally and axially excited states as discrete and continuous 
ones, respectively. As the typical temperature for the 
ultracold gas is on the order of the radial level splitting, 
the density of atoms populated in the radial ground state 
is obtained by solving the TBA equation, while those 
in the radially excited states are treated discretely as 
an independent ideal ID Bose gas in thermal equilib- 
rium. The density distributions of bosons and fermions 
are thus calculated in various experimental situations and 
the effects due to the interaction strength, particle num- 
bers and temperature are discussed. Our results show 
that the phase separation between bosons and fermions 
takes place at even lower temperature attainable by re- 
cent cooling techniques. 
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